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Abstract 

Determinantal point processes on a measure space (X, £, /j,) whose kernels represent 
trace class Hermitian operators on L 2 (X) are associated to "quasifree" density operators 
on the Fock space over L 2 (X). 

1 Introduction 

This contribution has been informed and inspired by several surveys of the topic of de- 
terminantal point processes that have appeared in recent years. 01 The first of these, 
Soshnikov (2000), is inspired by the determinantal point processes that arise in random 
matrix theory: the set of eigenvalues of a random matrix is a realization of a determi- 
nantal point process, if the random matrix is sampled from any of the unitary-invariant 
ensembles of Hermitian matrices (e.g., GUE), or from uniform measure on the classical 
(orthogonal, unitary, or symplectic) matrix groups, or from the Ginibre Ensemble. The 
review by Lyons (2003) is inspired by the Transfer Current Theorem jlj, which implies 
that the edges occurring in a randomly (uniformly) sampled spanning tree of a given 
finite graph G are a determinantal random subset of the edge set of G. Lyons's review 
concentrates on random subsets of countable sets, while Soshnikov's review is oriented to 
treat discrete subsets of a continuum. A very recent survey of determinantal processes 
(Hough et al. (2005)) includes the following newly- found example: the zero set of a power 
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series with i.i.d. gaussian coefficients is a determinantal point process [3] (the radius of 
convergence equals 1 almost surely). 

Hough et al. (2005) explain how a simple insight gives one a handle on number fluctu- 
ations in determinantal point processes. 0003111] The insight is that, in a determinantal 
point process with finite expected number of points, the distribution of the number of 
points is equal to the distribution of the sum of independent Bernoulli(Aj) random vari- 
ables, where < Xj < 1 are the nonzero eigenvalues of the "kernel" of the determinantal 
process. For example, consider the number of eigenvalues of a random n x n unitary 
matrix that lie in a given arc A of the unit circle. Denote this number by # n ^4- If the 
length of A is positive but less than 2n, then 

# n A - E# n A 

is asymptotically normal with unit variance. [71 1 101 ITT] The subset of eigenvalues that lie 
in A forms a determinantal point process on A, for it is the restriction of a determinantal 
point process on the whole circle, hence # n A is distributed as a sum of independent 
Bernoulli random variables. Thus, once one knows that the variance of j^ n A is (lnn)/7r 2 + 
o(n) \l'2l \1'6\ . the asymptotic normality of follows from the Lindeberg- Feller Central 
Limit Theorem. 

Determinantal point processes have a physical interpretation: they give the joint statis- 
tics of noninteracting fermions in a "quasifree" state. Indeed, this motivated the introduc- 
tion of the concept of determinantal (or "fermion") point processes in the first place. [H] 
Analogously defined "boson" point processes arise in physics and are called "permanental" 
point processes in probabilistic writing. |l4l !3] Recently, too, researchers have continued to 
investigate determinantal point fields from a quantum probabilistic point of view. [TBI I16j 
We adopt this viewpoint here, and realize that the satisfies of a determinantal point pro- 
cess with trace class Hermitian kernel JC on L 2 (X) are those of observables on the Fock 
space J-q(L 2 (X)) with respect to the density operator on J-o(L 2 (X)) that determines the 
gauge-invariant quasifree state with symbol K on the CAR subalgebra. However, we do 
not dwell below on the physical interpretation, nor do we discuss states on the CAR al- 
gebra in the following. Our main objective will be to construct the determinantal point 
process on X with kernel /C, when tC is the integral kernel of a Hermitian trace class 
operator on L 2 (X) with < < 1. Once the construction is understood, the fact that 
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the number of points in a measurable subset of X is distributed as a sum of independent 
Bernoulli random variables becomes obvious. 

Finally, let us remark that determinantal/permanental processes have a couple of 
different interesting generalizations. [T71 I18j And another rich survey of determinantal 
processes has just appeared in the electronic archive! [T§] 

2 Determinantal probability measures on finite 
sets 

Let X be a finite set, and let 2 X denote the set of all subsets of X. Let P denote a 
probability measure on 2 X , and let X be a random subset of X distributed as P. Then 
F(X D E) denotes the measure of the class of all subsets of X that contain the subset E. 
If there exists a complex- valued function K, on X x X such that 

F(X D {xi,x 2 , ■ ■ ■ ,x m }) = det (jC(xi,Xj)) ™ =1 (2) 

for all subsets {x\, ■ ■ ■ , x m } of X, where (JC(xi, Xj)}. . =1 denotes the mxm matrix 
whose (i,j) th entry is )C(xi,Xj), then P is said to be a determinantal probability 
measure with kernel K. The probabilities © determine the probabilities F(E) by 
inclusion-exclusion, hence there can be at most one determinantal probability measure 
with a given kernel /C. A very basic example of a determinantal probability on 2 X is the 
law of the random set produced by independent Bernoulli trials for the membership of 
each element of X; in this case the kernel K(x',x) = 5 x ' x F(x £ X). 

Suppose P is determinantal with kernel fC. Then the complementary probability mea- 
sure 

p c (X = S) = F(X = X\S) 
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is determinantal with kernel 2 — fC, where I(x', x) = 6 x 'x- To prove this, use the identity 

det (l(xi,Xj) - /C(x i ,x j )) i J =1 

m 

= 1 " + Yj det ( /C ( x ^' x ii>))a,fee{i,2} 

j=l l<ji<32<m 

det (fC(xj a , ^3b)) a ,be{v ,2,3} 
+ ... + (-irdet(/C(x,,^))J =1 . (3) 

The determinants on the right-hand side of (jSJ) are probabilities according to © , therefore 
det (l(xi,Xj) - K{xi,Xj))™: =l 

m 

= l-^PfeJcI)! E P({x il5 x i2 }cX) 

j=l i<jl<j2<m 

+ ... + (-l) m P({x 1 ,...,x m }cX) 
= P(X C <f \ {xi, . . . , x m }) [by inclusion-exclusion] 

= P((*\X) D{x!,...,x m }) 

= F(lDK...,x m }) . (4) 

Suppose that <Y = {xi, X2, • • • , x n } is an n-member set. Define the matrix Kij = 
(/C(xj, Xj)) i J_ 1 . If JT is a Hermitian matrix, then both K and I — K must be nonnegative 
matrices, since all of their submatrices have nonnegative determinants by © and 1)314(1 . 
Hence, if /C is the kernel of a determinantal random set and K is Hermitian, then if 
must be the matrix of a nonnegative contraction on C n , i.e., necessarily < \\K\\ < 1. 
Conversely, if K is the matrix of a nonnegative contraction on C n , then we will show that 
there exists a determinantal probability measure on 

2 {i,...,n} with kemel = Kij. 

The rest of this section is devoted to the construction of a determinantal probability 
measure whose kernel is a nonnegative contraction. Our point of view is that there exists 
a density operator on the Fock space over C n whose diagonal elements in the standard 
Fock basis give the desired probabilities. 

A density operator is a nonnegative Hermitian operator of trace 1. 
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Let J-(C n ) denote the exterior algebra over C n , i.e., 

F(C n ) = C © C™ © A 2 C n © • • • © A n_1 C n © A n C n , (5) 

where A m C n denotes the m th exterior power of <C n . The exterior algebra J-(C n ) is spanned 
by vectors of the form v\ A v 2 A . . . A v m , where V\, . . . , v m are any m vectors in C n and m 
is any number between 1 and n (together with an extra "vacuum vector" f2 to span the 
first summand). The expression v\ A v 2 A . . . A v m for vectors is formally multilinear in 
vi, . . . ,v m and satisfies 

Vj A • • • A Vi A • • • A v m = — v% A • • • A Vj A • • • A v m 

for j = 2, . . . , n. The exterior algebra J-(C n ) is 2™ dimensional and supports the inner 
product 

(vi A ••• Av m i, wi A • • • A w m ) = 5 m : m &et{{vi,Wj))™ =l 

(the vacuum vector is orthogonal to all v% A ... A v m and has unit norm). It can be shown 
that J-(C n ) is isomorphic to a subspace of the Fock space 

T (C n ) = CffiC ri ffi(C n ©C n ) ©•••©(©"£:") 

via the map that assigns 1 © 0c« © • • • © O^nc™ to Q and 

Oc © • • • © Og|m-l C ™ © S£[VI, . . . , V m ] © O m+l C n © • • • © 0®n C n (6) 

to v\ A V2 A . . . A v m . In ©, S£[vi, . . . , v m ] denotes the Slater determinant 

S£[vi,...,v n ] = —= S~] sgn(?r) U^(vi <g>v 2 ® ■ • ■ ®v m ) , (7) 

vreo m 

where S m denotes the group of permutations of {1, ... , m} and U n is the unitary operator 
defined on © m C n when ir E S m by the condition that 

U w (wi © W 2 © • • • © W m ) = ^-i(i) ® ^-1(2) ® • • • 8> ^Tr-l(m) ( 8 ) 

for all w\, . . . ,w m S C n . Henceforth, we identify the exterior algebra J-(C n ) with this 
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subspace of J-Q(C n ), and call it the "fermion Fock space." 

An orthonormal basis of J-(C n ), called a Fock basis or "occupation number" basis, 
can be built using any ordered orthonormal basis v = (vx, . . . , v n ) of C n . The vectors of 
the Fock basis can be conveniently indexed by subsets of {1, . . . , n}: the empty subset 
of {1, ... , re} corresponds to the vacuum vector Q and a nonempty subset . . . ,j m } C 
{1, ... ,n} with ji < ... < j m corresponds to the vector Vj t A ■ ■ • A Vj m . That is, the 
orthonormal set {f v (S) \ S C {1, . . . ,re}} is a basis for J r (C n ), where / v ({}) = and 
f v (S) = Vft A • • • A Vj m when S = {ji, . . . ,j m } with ji < . . . < j m . 

Suppose K is a nonnegative contraction on C n and let v = (vx, . . . , v n ) be an ordered 
orthonormal basis of C n such that Kvj = XjVj for all j. Let Dx denote the density 
operator 

D * = E { II " A ^)} (fyW' ■ )fy( s ) ( 9 ) 

5c{l,...,n} keS k£S 

on JT(C n ), where (f v (S), ■ )f v (S) denotes the rank-one orthogonal projector onto the 
span of f v (S). 

Proposition 1 Let K be a nonnegative contraction on C™ and let Dk denote the asso- 
ciated density operator Jj^ on the Fock space F(C n ). Then, for all ordered orthonormal 
bases w = (wi, . . . , w n ) of C n , 

S (f w (S),D K f w (S)) (10) 
is a determinantal probability measure on 2^ 1 ' , "' n ^ with kernel lC(i, j) = (Kwi,Wj). 

Proof: We first define the "second quantization" maps from operators A on (g) m C n 
to operators T m L4] on the Fock space J-o(C n ), and the dual maps from density operators 
D on J r o(C n ) to rre-particle "correlation operators" K TO [D] on ® n C n . 

Let i7(m, k) denote the set of injections of {1, . . . , m} into {1, . . . , k}. The cardinality 
of J^(m,k) is k^ = k(k — l) ■ ■ ■ (k — m + 1), the m th factorial power of k. For any operator 
A on (g> m C n , and any injection j £ J{vn, k) with k > m, we define the operator 

Ai3} = ^(lii)(2j 2 )---(rnj m )(^^ I ®---® / ) C/ (lii)(2i2)---(rnj m ) 

on (g> fc C n , where ^ r (iji)(2j 2 )-- (mj m ) denotes the permutation operator (JSJ) for the product of 
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disjoint transpositions (lji)(2j2) • • • (mj m ). Define r m L4] on JF (C n ) by 

r m [A] = o c e-90 r -i c »e Y A ij) ®---® Y Aii) ■ 

j£j(m,m) j£j(m,n) 

A density operator D on .F(C n ) extends to a density operator D® on J-o(C n ), which we 
will denote by D as well. The map A i — ► Tr(Z?r m [A|) is a linear functional on the space 
of linear operators on (g) m C n . Therefore there exists a unique operator K m [D] on ® m C n 
such that 

Tr{T m [A]D) = Tv(AK m [D]) 

for all linear operators A on ® m C n . In physics language, K m [.D] is the m-particle corre- 
lation operator for the state with density operator D. If Dk is defined as in Q, the key 
identity 

m times 



K m [D K ] = K®K®---®K sgn(vr)^ (11) 

may verified by comparing matrix elements of both sides with respect to the basis {vj l ® 
■"®Vj m \ 31,- -,3m G {1, ...,n}}. 

Given an ordered orthonormal basis w = (w\, . . . , w n ), let PJ denote the projector 
{wj, -}wj for j = 1, . . . , n. For distinct x±, . . . , x m , the operator T m [P™ (g) • • • ® ] is 
diagonal in the Fock basis {/ w (-5) I 5 1 C {1, . . . , n}}, and 



r m [p x w 1 ®...®p-j/ w ( < s) 



/w(S') if {x\,..., x m } c 5 
otherwise. 
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Let F(X = S) denote the probability (JTUJ). Then 



P(X D {x 1 ,x 2 ,...,x m }) 



SD{xi,X2,...,x m } 



£ (U(S),D K f w (S)) 



S^{xi,X2,...,x m } 



Tr((P™®---®P:jK m [D K }) 



det ({Kw Xi ,w x 



det (/C(. 




This proves the proposition. 



□ 



3 Determinantal finite point processes 



A finite point process on is a random finite subset of a space Let £ be a cr-field of 
measurable subsets of X.2Q A finite point process on (X, S) is specified by the probabili- 
ties Po,pi,P2, ■ ■ ■ that there are 0, 1, 2, . . . points in the configuration, and, for each n such 
that p n 0, a symmetrical conditional probability measure p n on (X n , Now let 

\x be any positive "reference" measure on (X, £). A finite point process is determinantal 
on (X, S, /i) with kernel K : X X X — ► C if 



for all disjoint, measurable Ex, ... , E m , m > 1. [3] If /C(x, y) is the standard version of the 
integral kernel of a nonnegative trace class contraction K on L 2 (X, S, /i), then there exists 
a unique [21] determinantal point process on (X,T,,fj,) with kernel IC(x,y). Conversely, 
if the kernel of a determinantal point process on X is the integral kernel of a trace class 
Hermitian operator K on L 2 (X), then must be a nonnegative contraction. |23| 

In this section we construct the determinantal point process on X whose kernel is the 
standard kernel (|17|) of a given trace class operator K on L 2 (X) with < ||ET|| < 1. This 
is accomplished by constructing a density operator on ^ r (L 2 (Af)) as we have done in the 



E(n#(*ni3,-)) = / ••• / det(K(x i ,x j ))™ =1 »(dx 1 )---Li(dx m ) (12) 
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preceding section — our quantum probabilistic point of view. There are many other ways 
to accomplish the same end, with or without our point of view. The original approach of 
Macchi (1975) was to start with a formula for the Janossy densities |24| of the desired point 
process, and then to verify (|12|) for that process. Soshnikov (2000) attacks the problem 
by first showing that certain Fredholm determinants involving K define factorial moment 
generating functions for finite families of random variables {^{X n Ej)\Ej E £}, and 
then constructing the desired determinantal point process via Kolmogorov extension from 
its finite dimensional distributions. Lyons (2003) uses the geometry of Fock space, but 
only in the case where K is a finite rank projector, then dilates nonnegative contractions 
to projections on a larger space to handle the general case. Hough el al. (2005) verify 
directly that kernels of finite rank projectors yield determinantal point processes, then 
treat the general mixture, in the probabilistic sense, of determinantal processes 

with projector kernels. 

Any density operator on T{L 2 {X, S, fi)) of the form D = ®D n defines a finite point 
process on X as follows. p n = Ti(D n ) is the probability of the event the configuration has 
exactly n points. The measure p n is absolutely continuous with respect to ® n /i and it is 
defined by way of the isomorphism 

® n L 2 (X,Y,,p) = L 2 (;t"\®"£,®». 

Regarding D n as an operator on L 2 (X n , (g) n E, <8> n /u), define p n {E) = p~ lr Tr(D n M.E), where 
Me denotes the operator on L 2 (X n ) of multiplication by the indicator function of E E 

Given a trace class nonnegative contraction K on L 2 (X), a density operator Dk on 
J-(L 2 (X)) may be defined using the spectral information in K as was done in Q above: 

oo 

d * = e e { n a * - a *)} ins), ■ )f(s) , (13) 

n=0sc{i,2,...}: fceS k&S 
#S=n 

where Ai > A2 > • • • are the eigenvalues of K and {f(S)} is the Fock basis constructed 
from the eigenvectors of K (really, any extension of an orthonormal system of eigenvectors 
of K to an orthonormal basis of L 2 (X)). Dk has the form (B(Dk)u and it can be verified 
|25| that the m-particle correlation operator K m [-Dx] exists and satisfies the key identity 
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Let p n and p n be as defined above, for Dk- Let E denote the expectation with 
respect to random point process defined by these p n and p n . Then 

m 

e( [] #(X n E,-)) = Tr(r m [Px ® • • • ® P m p/r) , (14) 

3=1 

where Pj denotes the orthogonal projector on L 2 (X), for both sides of (jTl)) equal 

E J^n njli #({^l> • • • !^n} n E j )p n (dx 1 ■ ■ ■ dx n ). But 

n>m 

Tr(r m [Px ® • • • ® P m ]^) = Tr ((Pxif ® • • • <g> P m K) ^ sgn(Tr)C^) , (15) 

7r£5 m 

since K m [D^] = ( ® m K) ^ sgn(-7r)f7 7r . Finally, one may verify that 

7r£S m 

Tr(( ®™ =1 P d K) sgnW^) = / ••• / det(AC(z i ,a; i )) f ™ 1 //(da; 1 )---//(da; m ) 

(16) 

if /C(x, y) is the usual version of the integral kernel of K, i.e., 

£(x,y) = ^2Xj<f>j(x)<f>j(y) (17) 

3 

where K(j)j = \j(f>j and ]TJ Aj = TrK. Equations (|14l) - (|16|) imply Q12|) holds; the finite 
point process defined through Dk is determinantal with kernel /C. 

Now that we have constructed the process, we can see immediately from that the 
total number of points in a random configuration is distributed as the sum of Bernoulli (Xj) 
random variables. In particular, n(l — ^fe) * s the probability that there are no points at all. 
This equals the Fredholm determinant Det(J — K). Let E be a measurable subset of X . 
It is not difficult to check via (|12j) that the determinantal point process on (E, T,\e, (i\e) 
with kernel Ke = McKMc\ E is the restriction to E of the determinantal point process 
with kernel K on X. Hence the probability that there are no points in E equals the 
Fredholm determinant Det(J — Ke)- In the context of random matrix theory, this yields 
formulas for the spacing distributions of eigenvalues. |271 I2H] 
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In case \\K\\ < 1, set L = (I — K) 1 K. It is easy to check from that 
(D K )n = ^ Det(I-K){® n L} sgn(vr)^ 

by comparing matrix elements of both sides of this identity with respect to an eigenbasis 
of K. This identity yields the determinantal formulas for the Janossy densities. j!4[ I21| 

4 Determinantal processes of infinitely many points 

Suppose that X is a locally compact Hausdorff space satisfying the second axiom of 
countability [2J2|, and let E denote the Borel field of X. In this context, a point process is 
a random nonnegative integer- valued Radon measure (a Radon measure is a Borel measure 
which is finite on any compact set). [3] Let /ibea cr-finite Radon measure on X .[S] A point 
process on (X, is determinantal with kernel /C if Q12|) holds. 

Most work on determinantal processes with infinitely many points has been done for 
the cases where X is a countable set with the discrete topology and \x is counting measure, 
or X is a connected open subset of M rf and \x is Lebesgue measure, or X is a finite disjoint 
union or Cartesian product of said spaces. If K is a locally trace class Hermitian operator 
on L 2 (X) such that < \\K\\ < 1, then (a version of) its integral kernel is the kernel 
of a determinantal point process on X.\Q This point process is the limit in distribution 
of the determinantal processes with kernels lc(x))C(x, y)lc(y), where C ranges over an 
increasing family of compact subsets of X. Conversely, if the kernel of a locally trace 
class Hermitian operator K defines a determinantal point process, then K must be a 
nonnegative contraction. [Tl HI] 

The case of a countably infinite set X with counting measure is treated in detail in 
Lyons (2003). In this pleasant special case, Hermitian operators on L 2 (X,2 X ,#) are 
automatically locally trace class. On the other hand, equations © and ©-© readily 
imply that K and I — K are both nonnegative operators. Therefore, the kernel of a 
Hermitian operator K on L 2 (X, 2 X , #) is the kernel of a determinantal point process on 
(X , 2 X , #) if and only if K is a nonnegative contraction. 
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